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o _ 

I Abstract. We consider an integrable system of reduced Maxwell-Bloch equations that 

^ ^ ■ describes the evolution of an electromagnetic field in a two-level medium that is inhomoge- 

^ i neously broadened. We prove that the relevant Backlund transformation preserves the reality 

of the n-soliton potentials and establish their pole structure with respect to the broadening 
parameter. The natural phase space of the model is embedded in an infinite dimensional 
QQ ■ loop algebra. The dynamical equations of the model are associated to an infinite family 

of higher order Hamiltonian systems that are in involution. We present the Hamiltonian 
^ functions and the Poisson brackets between the extended potentials. 

OO I Key words: loop algebras; Backlund transformation; soliton solutions 

2000 Mathematics Subject Classification: 37K10; 37N20; 35A30; 35Q60; 78A60 

^ ' 

^ ' 1 Introduction 

Integrable systems are closely related to the inverse scattering method that serves as a means of 
' integrating the initial value problem, and also to infinite dimensional Lie algebras, else known as 

■ Kac-Moody Lie algebras or loop algebras. Studying the underlying loop algebra reveals encoded 

^ I properties of the equations that are inherited from the integrable character of the model. The 

Adler-Kostant-Symes (AKS) theorem gives the Lie algebraic formulation of the dynamics of 
the system. 

In [l], several nonlinear partial differential equations were realized as the compatibility con- 



I dition of two linear matrix systems of the form, 

-iX q{x,t)\^ ^ (A 

where q and r are the potentials satisfying a certain nonlinear evolution equation and A,B,C 
are functions of x and t. The Zakharov-Shabat dressing transformation [2] was then employed 
to produce the soliton solutions of the system. The matrices appearing on the right-hand-side 
of ([1]) are 2x2, traceless matrices that can be viewed as elements of the finite dimensional 
Lie algebra 5 [(2). 5 ((2) may be prolonged to an infinite dimensional Kac-Moody Lie algebra 
with the aid of the spectral parameter A, and lead to an infinite number of systems in invo- 
lution. This is made exact in the context of the AKS theorem. Briefly, the theorem states 
that if one starts with a set of commuting functions on a Lie algebra, then the corresponding 
Hamiltonian systems are of course trivial. However, if we project those functions to appropriate 
subalgebras then the resulting Hamiltonian systems need not be trivial and continue to be in 
involution. 

In [SlIllISlElITlIElin], several integrable equations have been studied in the context of the 
inverse scattering technique and the AKS theorem. For example, the Toda system and the 
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Korteweg-de Vries equation, serving as representatives of the ordinary and partial differential 
evolution equations respectively, have been associated to Lax pair equations with one degree 
of freedom represented by the spectral parameter. The phase space of the relevant model was 
extended to include infinite dimensional loop algebras. The loop algebra was then decomposed 
into a vector space direct sum of subalgebras, and with the aid of the trace functional it was iden- 
tified with its dual. In that way, the systems obtain the Kostant-Kirillov symplectic structure, 
and an application of the AKS theorem revealed an infinite number of integrable Hamiltonian 
systems in involution. 

In this paper we consider a reduction of the Maxwell-Bloch equations that models the op- 
tical pulse propagation of an electric field through a two-level medium in the presence of an 
external constant electric field. The optical resonance line of the medium is inhomogeneously 
broadened. Following the terminology used by McCall and Hahn [10] and Lamb [11], we refer 
to inhomogeneous broadening as the phenomenon that occurs when the atoms of the medium 
possess different resonant frequencies due to microscopic interactions between them. In such 
a case, the induced electric dipole polarization is represented as a continuum, and the resulting 
optical resonance line is inhomogeneously broadened. In solids, such a broadening could be 
caused by a distribution of static crystalline electric and magnetic fields and in gases by the 
distribution of Doppler frequencies. 

Since the late sixties and seventies with the papers [10 ^ 111 ^ 112]. among others, the Maxwell- 
Bloch equations have undergone several treatments. Recently, various reductions of the equa- 
tions have been studied both analytically and numerically. Lax pair operators, Darboux trans- 
formations and soliton solutions were constructed and analyzed [13l [HI [151 [E] and interesting 
applications in crystal acoustics have emerged |17[ 118] . 

Our scope in this paper is the study of the integrable structure of a reduced Maxwell-Bloch 
system and the connections that arise with Kac-Moody Lie algebras. In particular, we prove 
that the Backlund transformation preserves the reality of the n-soliton potentials Vn G N, and 
establish their pole structure with respect to the broadening parameter. The solitonic phase 
space of the model is embedded in an infinite dimensional loop algebra and an application of 
the AKS theorem allows us to view the system as a member of an infinite family of systems in 
involution. We present the higher order Hamiltonian functions and flows, as well as the Poisson 
brackets between the extended potentials. 



2 Phase space 

The optical equations we shall consider are the ones presented in [19] . They model the propa- 
gation of an electric field in a two-level quantized medium, where the optical resonance line of 
the medium has been inhomogeneously broadened. The classical wave equation of Maxwell ([2]) 
is used for the evolution of a unidirectional electric field and is coupled with the quantum 
mechanical Bloch equations ©-([SI), that describe the behavior of the induced polarization 
field, 

de de , ^ , 

^ + ^ = (^5.)^ , (2) 
^ = {f3-ie)S^, (3) 
^ = -(/3-7e)i?. + -a;C/^, (4) 
^ = -2uS^. (5) 
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{f)g = f{uj)g{uj)duj, and denotes the weighted average of the function /(w) with respect to 
the distribution function, 



7r((a; - Wq)^ + d"^) ' 

For a physical interpretation of the model see [19]. In this paper we shall study the system 
(HI)"® from a Lie algebraic point of view. 

The system is completely integrable and admits a Lax pair representation. We define the 
differential operators L and A whose commutativity, [L, A\ := LA — AL = 0, is equivalent to 
equations ([S])-® 

where, 

= X{hon + foT) + eo£, 

/oo 1 
r^[X{hin + fiJ^) + ei£]duj, 
-oo — ^ ) 

and 

ho = ^, hi = -^-fiog{u;)R^, 

fo = 0, /i = -^jujg{uj)S^, (6) 

eo = - 7e), ei = l-fuj^g{uj)U^. 

TC, J- and E form a basis of the semi-simple Lie algebra su(2), and are given as follows, 



We call /ij, fj, Cj for j = 0, 1, the potentials and Q^^\ Q^^'^ loop elements because they can be 
considered as elements of an infinite dimensional loop algebra that we will define in Section [3l 
We note that the potentials depend on the solutions e, R^^, S^j, U^^ of the inhomogeneously broad- 
ened reduced Maxwell-Bloch (ib-rMB) equations ([2])~©- The commutation of the differential 
operators L and A gives rise to the following Lax pair equation 

which is equivalent to the ib-rMB system. 

The Lax pair can be used to construct a Backlund transformation (BT) that iteratively 
produces the soliton solutions of equations ([2])-®. 

We consider the spectral problem, 

and aim to find a new eigenfunction ^' and the corresponding new loop element Q^^^ that satisfy 
the spectral problem. The loop elements are functions of hj, fj, ej, for j = 0,1 and will, in turn, 
give rise to the new solutions of the ib-rMB system via expressions ®. This transformation 
theory leads to an analogue of superposition formulas that allows one to construct multi-soliton 
solutions starting from single solitons by algebraic means \21\ [22| [23l [2^ [T5] . We briefly describe 
the procedure and quote the relevant theorem from reference |19j . 
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One begins with a constant solution to equations ©-([S]), which in turn determines poten- 
tials ^ and the corresponding loop element, call it Qq. We then find a simultaneous, fundamen- 
tal solution ^'i to the Lax pair system = 0, = and define = ( ) := ^'i(A = i'i)c{, 

where c{ = i ^\ ) is a constant vector with c\,c^ G M, and the matrix A^i as, 



iVi 




^2 



The BT matrix function is constructed as: 

Applying G to ^'i yields a new fundamental solution: ^'2(2^1, ci; A) = cl; A)^'i(z^i, cl; A), 

and the procedure is iterated. The formula for the loop element after n iterations of the BT, call 
it Qn, in terms of the previous one Qn-i and the matrix Nn is the context of the next theorem. 

Theorem 1. 



-( 



(cj2 - A2) (cu2 + ml) 
X {A [i^^hJl-^-H + fr'^) - {mn)\Nn'HN-^){Kl-^n + fr^F)[Nn'HN-^) 

+ Loh^-^S - {mnferHNnnN~')£{NnnN-')} duj. (8) 

We have taken the specific value of the spectral parameter to be purely imaginary, = 
irrin G iM. The general form of the n-soliton loop is given by, 



Qn{\) = XihoH + /o.F) + e'^S + r — 1— -[A(/i?H + f^T) + e'l8]dLo, 

J -00 (w^ - A2) 



(9) 



where the upper index of Cq , /i" , /f , e" and the lower index of Qn, ^n, ^n, indicate the level 
of the Backlund transform. We note that /iq and /o are constant functions of space and time and 
are invariants of the level of the BT. To ensure the reality of the potentials and consequently of 
the solutions e, R^j, S^j, U^j, we compare ([8]) and ([9]) and impose the following conditions: 

a) [n, NnnN-\ [T, NnnN-\ {NnnN-^)£{NnnN~') € span {£}, (10) 

b) {NnnN~^)n{NnnN-^), {NnnN-^)T{NnnN-^), [e, NnUN;;^] e span{H,^}. (ii) 

By definition 

Nn = Re{(t)i )T + lm{(t/l)n - Re((/.^)^ + Im((/)^)5, 

where Z is the 2x2 identity matrix. One can find that conditions (|10 p - (|lip are satisfied if and 
only if 

Im(,/.^)Im(</.^) + Re(0^)Re(,/.^) = 0. (12) 

We can construct (j)^, cj)2 such that one of the following two cases holds: a) G M, 02 G iM, or 
b) G iM, 4>2 £ I^- This is the context of the following proposition. 
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Proposition 1. At any given level of the Bdcklund transformation the reality of the n-soliton 
potentials can he secured by an appropriate choice of the transformation data. 

Proof. A constant set of solutions to equations ©-(in]) is given by e = ^, 5^ = 0, C/(^ = 
and = a nonzero constant. The corresponding potentials become Hq = ^, hi = 

— ■^'yujg{u!)R™^^, /o = /i = eo = = 0. Following the procedure of the Backlund transform we 
find, (pl = Re(c];)e^'^ +ilm{c\)e^^ , 4)\ = ^e{c\)e~^^ + ilm{c\)e~^^ where xi = 2{{vi + m.iho)C — 
mih^r) G M, and vi is an expression independent of uj, A, C, and r. We choose cj € M and c\ G lU. 
so that (/>J G M, (/>2 € iM and condition ()12p holds for n = 1. The proof that the condition is 
satisfied at any given level of the BT lies on the following: 

N2n-i € span{X, £■}, and G span{H,!F}. This yields Gn(A) G span{A/, Vn € N. 



Using the definition = G'„(A = im„)4'„_i(A = we deduce that ^'2n-i S 



and ^ I [^]^ ^jgj J- Therefore, if we choose c„ € <j ( ) and use = ^nCn, we obtain 
that ^2n-i £ I ( ^ ) r 5 s-iid $2n S \ (^mi \ desired 



We aim to give an integral- free representation of Qn{^) so that its A-structure becomes 
apparent. We begin with Qi(A) and then generalize the construction for Qn{}^)- The general 
form of Qi{X) is the following, 

/OO 1 
-2- [Kh\H + flJ') + e\£] dio. 
-OO {^^ - A^j L 

The one soliton potentials are given as, 

1 1 2miLo'^hi 

en = 2mi/io sech(xi), Ci = ^sech(xi), 

(w^ + mf) 

/i = / 9 ^ — 9T sech(xi) tanh(xi), /i} = —5 ^ Tu;^ +m{- 2m\ sech^(xi)l , 

(cj^ + my (w^ + m|j 

where xi = 2{{vi+miho)Q-mihQT)+\xi{c\/cl), hi = -\^ujg{uj)R-^'' and vi = -^ (g'+aFCg) - 
We write the last three in the following more convenient form, 

1^ Pl^ ..ON 

(u;2+^2)((^_^^)2+^2)' ^ 

fi = n4i 

(u;2 + m2)((a;-a;o)2 + a2)' ^ ^ 

1 (a;2 + „if)((a;-a;o)2 + a2)' ^'^^ 

where /32, /^a, /?4 are functions of and r, but do not depend on uj or A. We will use 
contour integration and Cauchy's integral formula to compute Qi. We define the complex- 
valued function, 

h{z) = 2 I ^^{A[(/?3Z2 + /34)W + /32^] + 

(z^ — A^)(z^ -I- m|)((z — Wo) + c"" j 

which has three poles in the upper half complex plane at zi = X, Z2 = imi, and z^ = loq + ia. 
We integrate around a simple contour that consists of a semicircle of radius R in the upper-half 
complex plane, call it Cr, and the segment on the real axis from —R to R. We choose R big 
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enough to include all the poles of h{z) that appear in the upper-half complex plane. It is not 
hard to see that lim h{z)dz = 0, and thus using the Cauchy integral formula we obtain, 

3 

Qi = A(/ioW + /o^) + el£ + 27ri ^ Res^=,,/i(z), 

k=\ 

which can be written as, 

Qi(A) = \{h^H + foT) + el£ + —^-^[X{6in + 82^) + 638] 

(A^ + mf) 

x{6\n + (5» + + xH5\n + Is-^) + x^s^s 

where 6j, 5j, 6j, j = 1, 2, 3 are independent of A and lo. Let Si = {61,62,63}, S2 = {6162, 63} and 

S3 = {6162, 63}. The elements of each set Si, S2 or 5*3 are functions of the one-soliton potentials 
and they obey the ib-rMB equations. 

We will use induction to obtain the general form of the n-soliton potentials /i", e". 

Proposition 2. The pole structure with respect to the broadening parameter lo of the general 
n-soliton potentials of the ib-rMB system is the following: 

n 

Ln K=u fn _ k=0 

"1 — n ' Jl 



n 

k=0 


((a; 


n 




k=l 




n 

uj E efct^^'' 




k=0 


{{lo 


n 

-a;o)2 + cT2) \[{LO^ + ml) 

k=l 



{{u;-LOoy + a^)Ui^' + rnl) 

k=l 



We note that 'yk,6k,ek, k = 1, . . . ,n, are analytic functions of C and r, and do not depend on uo. 

Proof. The proposition holds for n = 1 as can be seen by ()13p -()15 p . We assume that the 
proposition holds for n — 1 and show that it holds for n. Using the induction hypothesis in 
Theorem [1] and the fact that NkHN^^ G span{7Y, J^} for any /c G N we get, 



Qn{\) = \{hon + foT) +e^£ + j 



CO 



^ (w2 _ ^2)((^ _ ^q)2 + ^2) Yl (c^2 + ^2 



k=l 

( r "n. "In "\ 

X i A " ^ 



.k=0 k=0 



+ ^ ekLo'^'^S > dLo. 



k=0 



By definition 

Qn{\) = X{hon + foT) + e^^ + —L^[x{Hin + /i"^) + e^^£]duo. (17) 

J-00 {LO — A ) 

Therefore, by equating the last two expressions of Qn{X) we prove the proposition. We note 
that some of the coefficients 7^, 6k, Cfc may be equal to zero. For example, for n = 1, 61 = 
and eo = 0. ■ 
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The form of the n-soliton potentials can be used to identify the A-structure of a general n-soliton 
loop element Qn(A). 

Proposition 3. 



Qn{\) = KhoH + foT) + e^£ + 



((u;2 - a2 - A2)2 + 4w2^2 



+ E (A2 + ml) [^^'^'^ + '^'•^^ + ^^^^ 

The proof follows the same idea that was used to derive the integral-free form (jl6p of Qi (A) and 
is omitted. 

Let Sk = {6f,62,S^}, k = 1, . . . ,n, Sn+i = h} and Sn+2 = h}- The elements 

of each set 5^, k = 1, . . . ,n + 2 are functions of the m-soliton potentials for m = 1, . . . ,k and 
they satisfy the ib-rMB equations. 



3 Embedding 

We call Ok = irrik, k = 1, . . . ,n, a„+i = loq — ia and 0^+2 = "^o + ^f- Then Qn{^) of (fTH]) can 
be rewritten as 

n+2 + 



Q„(A) = A(/io7^ + /o^) + eo^ + V + 

f— ^ A - Q!fc A + afc 
fe=i 



where , are linear combinations of the basis elements of the Lie algebra su(2), {Ti, £}. 

We define the following infinite dimensional Lie algebra Q^^^ that extends loops of the 
form (fTHI). 



{00 n+2 1 



We note that the coefficients hj , fj , of ([T9|) are not the same as those of pT|) where the upper 
index denotes the level of the BT. In (1191) X is an element that extends the natural solitonic 



phase space of the ib-rMB equations and the upper index indicates the relevant pole ak, whereas 

4- 



the lower index indicates the order of the pole 



We embed this Lie algebra into a larger one, 

Q-- 



2(n+2) 



where we have set ak+n+2 = — a/ci for A; = 1, . . . , n + 2. We note that for Q to be a Lie algebra 
a finiteness condition needs to imposed. Namely, Xj = and Xj = 0, V A; = 1, . . . , 2(n + 2) and 
j > jo for some jo € N. 



4 Application of the Adler— Kostant— Symes theorem 

Following the ideas in the theorem of Adler, Kostant and Symes [6l [71 [8] , we decompose the 
infinite dimensional loop algebra Q into a direct sum of two subalgebras a and b, and define an 
ad-invariant, non-degenerate inner product on Q. The perpendicular complements a"*", 6"*" with 
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respect to the inner product serve as another direct sum decomposition of Q. Using the Riesz 
representation theorem and the non-degenerate inner product one may define an isomorphism 
between and 6*, the dual of the Lie subalgebra h. The canonical Lie-Poisson bracket that 
exists on b* is then represented on . If necessary, one may translate a-*- by any element (3 £ b-^ 
that satisfies (/3, = \/ x,y £ a, so that the translated space P + a-^ includes the natural 

phase space of the relevant system [20] . 

We write Q as the vector space direct sum of the Lie subalgebras 

{oo 2{n+2) X'' ] f oo " 

^^^^ ^ = E E (A^tyj' b=}xeQ: X = E^^^' 

A non-degenerate inner product is defined on Q using the trace map. Namely, 

Q X Q ^ c, (X, y) ^ Tr(Xy)o, 

where {XY)q denotes the matrix coefficient of A*^ in the product XY. We note that the inner 
product is ad-invariant. That is, (X, adyZ) + (adyX, Z) = 0, where adxY = [X,Y]. The 
perpendicular complements of the subalgebras a and b with respect to this inner product take 
the form. 



oo 2(n+2) -^k 

XgQ:X = Xo + EE ' 



X eQ: X = ^XjX^ 



(A - OkV 
j=i k=i ^ ^' 

The natural phase space of the ib-rMB equations, as can be seen in (jlSp . contains elements that 
belong in as well as terms of order A^. Therefore we translate the space hy (3 = XXi € b^. 
We note that [X, y]) = VX, y G a, which is a necessary condition for the AKS theorem. 
We consider the set. 



+ /? = < 



oo 2(n+2) 

XeQ: X = X, + XX, + Y^^ Y: (A^tlV^ ^^l) 



which includes the phase space of the ib-rMB equations, and define a Lie-Poisson bracket on 
a"*" + f3. If $ and are functions on a-*- + /3 we first compute their gradients V<I>, V\I' in the 
full Lie algebra Q and then define 



{<^,^}iX) = -{X, 



llv^x),llv^ix) 



yXea^ + p. (22) 



We note that denotes projection on the Lie subalgebra b. VH{X) denotes the gradient of 
a function H on Q and is defined as follows, 

(y,VFm) = li,n^S^±in-S^). (23) 

We quote the AKS theorem and show how it can be applied in the case of the ib-rMB equations. 

Theorem 2. If H, F are invariant functions on Q* = Q, and we denote by H^^°\ F^^°^ their 
projection to the subspace a-^ + f3 then {i^P^j^ _pp™j} = Q, and the Hamiltonian system associated 
with the invariant function H is given as, 

^(/3 + a) = -ad^^VH(/3+a) + V/3 + ae/3 + a^. (24) 
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Remark 1. ad* is negative the dual of the ad map. 



Remark 2. The definition of H being an invariant function is ad^j^f^^^x = 0, Vx € Q. 

Remark 3. The existence of an invariant inner product allows one to identify the ad-map 
with ad*, because on one hand by definition (adxy,z) = —{y,ad*z), and on the other a simple 
calculation using the definition of the inner product shows that (ad^y, z) = —{y, adxz). 

If a set of invariant functions is given, the theorem guarantees their commutativity in the 
canonical Lie-Poisson bracket and gives the form of the Hamiltonian system associated with the 
invariant function. To apply the theorem we construct a set of invariant functions {<I>fe}fegN via 
the following operator, 

2(71+2) 

Mk{x) = n - "^)'^- 

i=l 

The invariant functions are defined as ^^{X) = ^{Mk{X), X). We compute the gradient of the 
By definition, 



(y,V$fc(X)) = lim 



•^kiX + eY) - ^k{X) 



^2(n+2) 

n {x-m)''x,Y 

i=l 



(y,Mfc(x)). 



Thus V$fc(X) = Mfc(X). According to Remarks 2 and 3, to actually demonstrate the invariance 

2(n+2) 

of ^>fc we must show that [V$fc(X),X] =0. We have that [V^>fc(X), X] = J] {X-ai)^[X, X]=0. 

i=l 

Therefore the functions <5fc, A; € N are invariant and can be used in the context of the AKS 
theorem, which reads as follows. 



dtk 



ll\7^k{X),X 



since = [V<Pk{X),X] = [Ua'^^kiX), X] + [l\b^<^k{X), X]. We used the subscript k for 
the time variable tk to distinguish between the different dynamical evolutions of the systems 
associated with the Hamiltonian functions For each G N the following systems are in 
involution: 



dX 
dtk 



llMkiX),X 



For A; = and a truncated X of the form (119^ where Kj, fj, e*- 



n+2 



X = XihoH + hT) +eo8 + Y,- 



the Hamiltonian system ()25p becomes 



a: 



(25) 

for j > 2 and i = 1, . . . ,?2 + 2, 
[\{h\n + f{T)+e^^8l (26) 



dt(\ 



(27) 



which is the zero curvature representation of the Lax pair equation ([7]) for the ib-rMB equations. 
We have thus identified the ib-rMB system as a member of the infinite family of systems 
that commute with respect to the canonical Lie-Poisson bracket (|22p . 
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5 Extended flow 

We consider a general element of the Lie algebra Q of the form, 

oo n+2 ^ 



or equivalently, 



j{n+l) 



n+2 

n - 



(29) 



We remark that /i*-, /j, e*- in (j28p are not the same as the ones in (j29p . The latter ones are linear 
combinations of the former. However, to avoid introducing yet another symbol we use hj"^'^^ , 
j2m.+i ^ g2m -j^ ([29]), where j indicates the order of the pole at = of and 2m + 1 or 2m the 
power of A multiplying Ti, T or £ respectively in the numerator. 

To examine the extended flow associated with the ib-rMB equations we set A; = in the 
involutive systems (j25p . The relevant system takes the form, 

dX 
dto 



where V^q{X) = Mo{X) = X. We write X as follows: 



.=0 n (A2 - ajy 

i=l 



where Pj is a polynomial in A given as 

j{n+l) 

Pj{X)= >?'^^\hf+^H + f]'^+^T) + \^'^ef^£, j = l,2,..., 

and Po = A(/io'H + f^T) + eof . The degree of Pj{X) is given by, deg(Pj(A)) = 2j(n + 1) + 1. 
Projecting V<^o(^) to the subalgebra h is equivalent to keeping the polynomial part of X, which 
is Pq. Thus the Hamiltonian flow for = takes the form, 



P 



' n+2 

n (A2 - o^iy 

i=l 



which unravels to 



dX 



oo nin+l) 



/'n+2 



j=l y m=0 J \i=l J 

The coefficients of the matrices 7^, and £" appearing in (j30p are defined as follows. 



(30) 
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de^,. = 2(/io/J"^+^-/o/if+^). 

We observe that the expressions multiplying Ti, and J- in (j30p have no polynomial part since 
the degree of the numerator that equals 2j[n + 1) + 1 is strictly smaller that the degree of the 
denominator that equals 2j{n + 2), for j > 1. However, the expression multiplying £ carries 
a polynomial term. In particular the evolution equation (jSOp can be written as 



dt 



,y oo /j(n+l) i(n+l) 

j=l \ m=0 m=l 

3{n+l) \ /n+2 \ 

+ E Cr^~i,~i)in+i)>?'^d4]^'h / n(^' - ^"y 

m=0-l)(n+l) / \i=l / 

oo /n+2 \ 

j=2 \i=l J 

n+2 

where Cn+i-m = (—1)'" ' ' ' o^i^+i ) ^ {l,...,n + 2}. On the other hand by the 
definition of X we have that, 

? = A(^W + 1°^) + 
ato "io "^0 dtQ 



3=1 \ / \*=1 / 

By equating the different powers of A we disclose the system induced by the Hamiltonian func- 
tion $0- As previously noted, for a fixed j G N, the elements of the sets Sm = {h?^'^^ , /J™""*"^ , e^™} 
satisfy the same system of equations for any m = 0, 1, . . . , j(n + 1). Therefore it suffices to con- 
sider only one such set. Without loss of generality we choose Sj(^n+i)- The evolution equations 
take the form: 

^ = 1° = 0, p= 2{hof!-+' - hh^r'). 

dto dto dto ^ ^ jv I J, 

' = 2(/oef - eo/f ("+^)+^), 



p2j(n+l)+l 



dto 

+ 2{Kf]^,^'^^^^-''^-'' - fohfl^^^^^^ 

To reveal the extended flow for the ib-rMB equations we set j=l: 

^ = 1^ = 0, P = 2(/io/^+=' - fohr^'), 
dto dto dto 

^t,2n+Z jr2n+3 

= 2(/oer+2 _ eo/2"+3)^ ^ ^ 2(eo/.r=^ " hoe'r'), (31) 
atfi ctrn 
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J 2n+2 /n+2 \ 



+ 2{hoft^'-foht^+'). 

We note that the couphng of the above system to the evolution equations satisfied by the higher 
order potentials that correspond to j > 2, is captured in the dynamical equation for e^""*"^. 
If /ig"^^ = _ then the system reduces to the dynamical equations that the n-soliton 

potentials of the ib-rMB equations satisfy. 

6 Hamiltonian functions and Poisson brackets 

In this section we aim to write the extended flow of the ib-rMB equations given by system (j3ip 
in Section [5l in the canonical Poisson form, 

^ = {eo,^o}, ^^ = {/^r^<fo}, (32) 

a /•2n+3 /?p2n+2 

^/l _ / f 2n+3 ^ X "^^1 _ r^2n+2 ^ i 

The Hamiltonian functions for the systems (j25p described in the context of the AKS theorem 
in Section m are defined as 

1 \ 

= -{M,{X),X) = -Tr n (A - c.^)'x'j . 

We let /c = and consider a general X of the form, 
X = \{hon + hF) + ee£ 



E 



OJ J J 



\ 



ri+2 
i=l 



The Hamiltonian function is found to be: 

^,{X) = \tt{X^)^ = el + 2{hohl^+' + /o/^+=^). (33) 
We define the following functionals, relevant to the potentials that appear in the Hamiltonian: 

\2n+3 

ho{X) = ho coefficient of XH, /if +^(X) = hf^^^ coefficient of — H, 

n (A2 - aj) 

i=l 

\2n+3 

/o(X) = /o coefficient of A^, fl^^^^X) = coefficient of — 

n (A^ - «D 

i=l 

\2n+2 

eo(X) = eo coefficient of e?"+2(X) = e?''+^ coefficient of — £. 

n (A2 - a?) 

i=l 



Prolongation Loop Algebras for a Solitonic System of Equations 



13 



Using definition (p3]) we compute the gradients of these functionals: 

1 \2n+3 1 

v^oW = -^^^ vh'r^^x) = --\n, 



a: 



2^ 



i=l 

1 \2n+3 1 



V/o(^) = J", V/^+3(^) = _1XT, (34) 

n - «?) 

1=1 

Veo(X) = -^£, Vef = -^A^^. 

For example, to obtain S/ho{X) we consider the equation 

e— >o e 

which implies that Tr(y V/io(X))o = ^o(^)- Therefore V/io(-^) is the element of Q such that the 
constant term (with respect to A) of the product YS/ho{X) has trace that equals precisely hQ{Y). 
We note that the matrices TiJ-, J-£, STC are traceless whereas Ti.^ = J-'^ = S'^ = —2. Having 
that in mind, we find that Vho{X) = — ^774:^ H. 

n (A2-a2) 

1=1 

The Poisson brackets between the potentials appearing in the Hamiltonian can be computed 
using definition (j22p . For instance, 



{eo,hl-+^}{X) = -lX, 



nveo(x),nv/if+3(x) 

b b 

1 . „\ „ / 1 



X, -\T j = -Tr = ^i"^'- 

In a similar fashion we obtain the rest of the Poisson brackets, 

/n+2 \ 
/n+2 \ 

Using the Poisson brackets (|35|) . we find that the extended flow for the ib-rMB equations given 
in (j3ip can be expressed as the canonical flow (j32p associated with the Hamiltonian function $o- 
Higher order functionals can also be defined using a diagonal formation that gradually sweeps 
all the potentials. In particular, for a general X of the form (j29p we define for N = 0, 1, 2, ... , 
the following higher order functionals (in bold to distinguish between ho, /q, cq and /iq, fo, cq): 

7V+1 N+1 



'J 

j=i i=i 

N+l 

2j{n+2)~2N~2 
j 



We note that if the upper index of the potentials that appear in the sums is less than zero 
then the potentials are set to zero. The set 5 = {eo, h]\f, ctv : = 0, 1, 2, . . .} includes all 
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the dynamical quantities that enter the AKS flows (j25p . The gradients of these higher order 
functionals can be computed using definition (|23|) . For instance, 

-, /ri+2 \ /r!.+2 \ 

1 /''^^ \ 



where V^o = V/l^"+^ V/q = V/l2"+^ Veo = Vef'^^, and are given in §4 

Working in a similar manner as in the example for {eo, /ii"^^}, one can find the Poisson 
brackets between the higher order functionals, i.e. {eo, hi}{X) = fi{X). 



7 Summary 

In this paper we have considered an integrable system of reduced Maxwell-Bloch equations, that 
is inhomogeneously broadened. We show that the relevant Backlund transformation preserves 
the reality of the n-soliton potentials Vn € N, and establish their pole structure with respect to 
the broadening parameter. We obtain a representation of the relevant phase space in the spectral 
parameter A, which is then embedded in a prolonged loop algebra. The equations satisfied by the 
n-soliton potentials are associated to an infinite family of higher order Hamiltonian involutive 
systems. We present the Hamiltonian functions of the higher order flows and the Poisson 
brackets between the extended potentials. 
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